Numerical simulation of Einstein-Podolsky-Rosen experiments in a local hidden 

variables model. 
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We simulate correlation measurements of entangled photons numerically. The model employed is 
strictly local. In our model correlations arise from a phase, connecting the electromagnetic fields 
of the two photons at their separate points of measurement. We sum up coincidences for each 
pair individually and model the operation of a polarizer beam splitter numerically. The results 
thus obtained differ substantially from the classical results. In addition, we analyze the effects of 
decoherence and non-ideal beam splitters. It is shown that under realistic experimental conditions 
the Bell inequalities are violated by more than 30 standard deviations. 
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Since experimental tests of the Einstein-Podolsky- 
Rosen (EPR) problem became feasible Q-^] , the field has 
gained a popularity and importance, which is startling. 
The experiments themselves do not seem to justify this 
broad interest. The fact that polarization measurements 
of two separated photons are connected, even if the mea- 
surements are in space-like separation ||, is in itself 
problematic, but hardly accounts for the large number 
of publications devoted to the problem every year. The 
importance of the EPR problem derives from its impact 
on two seemingly separate fields. On the one hand, the 
theoretical predictions impose a limit on theoretical con- 
cepts, aimed at modifying quantum mechanics by hidden 
variables models H0. Here, the current understanding 
is that the experiments contradict any local model ||Jj] . 
On the other hand, entanglement of separated photons is 
essential for many theoretical concepts in quantum com- 
putation and quantum information Jjoj . Here, the ques- 
tion whether quantum mechanics provides the only suit- 
able model for these experiments, is essential, because 
qubits, the principal information units in these concepts, 
are only defined within the quantum formalism. 

Theoretically, the entanglement of separated photons 
rests on the validity of the superposition principle. The 
question, whether superposition is a sufficient condition 
for non-locality, is open. Deutsch and Hayden showed 
recently that all information in an EPR experiment is 
stored locally, even though it cannot locally be accessed 
(□J. In a different approach it was shown that the ob- 
served connection can be interpreted as fourth-order op- 
tical correlation. The connection then arises due to the 
phase of the photons' electromagnetic fields (l^JT^]. All 
physical effects in this case are local, even though the cor- 
relations must be expressed in terms of polarizer settings 
at different locations. 

In this paper we present a model, where the phase 
connection between the two photons is retained, but the 
superposition principle for the separate photon fields is 
discarded. Essentially this means that interference terms 
in the correlations are omitted. We shall show, how- 



ever, that the model also in this case is in accordance 
with experimental results [0-0], provided the polarizer 
beam splitter (PBS) is not completely ideal. The non- 
ideal PBS in our model can only discriminate between 
the horizontal and vertical polarization of photons, if the 
intensities of the photon's electromagnetic field compo- 
nents, projected onto the orthogonal polarizer axes, dif- 
fer by more than 30%. We shall show that this threshold 
is consistent with and necessary for the high visibility 
observed in state-of-the-art experiments. Furthermore, 
we shall establish that the Bell inequalities are violated 
under these conditions by more than 30 standard devia- 
tions. Single photon experiments will be treated within 
the same model in a separate paper. 

The setup in our model consists of a source for entan- 
gled photons and two PBS with photon detectors. The 
source, at z = of the coordinate system, emits two 
entangled photons. Their planes of polarization are un- 
known. But it is known that these planes at the moment 
of emission are perpendicular to each other (AtyS = 7r/2). 
The polarizers are at z = ±L, where L is an integer mul- 
tiple of the photons' wavelength A. Since laser emissions 
have finite extensions in the range of a few microns [fli) , 
we describe each photon as an excitation of the electro- 
magnetic field of limited extension. To simplify the de- 
scription these fields shall be monochromatic. That is 
to say we disregard the effects of higher Fourier compo- 
nents, which arise from the limited extension of the fields. 
Then the phase along the photon's path of propagation 
is equal to the phase of a monochromatic wave. For pho- 
tons of circular polarization the angle of polarization at a 
point zq depends on the distance from the source and the 
initial angle of polarization. For linear polarization only 
the amplitude of the field depends on the distance Zq, but 
not the angle of polarization. This is changed, however, 
in EPR experiments using linear polarized light, by an 
optical modulator pq |. A modulator shifts one compo- 
nent of the electromagnetic field, either the horizontally 
or the vertically polarized one by a phase-angle Aip. If 
the field is initially in linear polarization, it is described 
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by E(z,t) — E cos(kz — u)t), where X — (E x,E 0Y ). 
The angle of polarization tpo = arctan {E v /E x ) in this 
case is independent of the phase. If one component, say 
E y is shifted by a phase-angle Aip, then the field after 
the modulator will be: 

E' = (Eox cos(fcz — ujt), E y cos(kz — ujt + Ai/>)) (1) 

The angle of polarization <p is then no longer constant, 
but depends on the phase: 

ip = arctan | 0Y [cos A-0 — tan(fa — tut) sin A-0] 1 (2) 
I E x J 

Therefore in both cases the angle of polarization ip de- 
pends on the location z of a measurement. The measure- 
ments are performed with polarizer beam splitters and 
two photon detectors JHJ. This differs from using opti- 
cal polarizers, because the combination PBS and detec- 
tor yields not a continuous result, e.g. the field intensity 
after the polarizer, but a discrete one. The photon is ei- 
ther detected at detector 1, or at detector 2. In one case, 
it is said to have vertical polarization (+), in the other 
horizontal polarization (-). The representation of this 
experiment by a classical optical model is not straight- 
forward. In Furry's integral, for example, the correlation 
is described by Jl6| , 

P(a,/3) = J dip cos 2 (<p + 7T/2- a) cos 2 (<p - (3) (3) 

where P(a, (3) is the probability of a coincidence, tp 
the unknown angle of polarization, and a and j3 are the 
polarizer settings. In this equation it is assumed that 
a single measurement can be described by the product 
cos 2 {tp + tt/2 — a) cos 2 (ip — (3). But this is evidently not 
the right representation for the dichotomic results ob- 
tained in the experiments. In our view it is partly this 
error, which is responsible for the deviation between the 
correlation curves obtained in the experiments, and the 
correlation curves obtained from Eq. (0). Analyzing the 
measurements event by event only two results are ob- 
tained in the experiments on either PBS: (+), if the angle 
of polarization ip is between 7r/4 and tt/2, and (-), if it is 
between and ir/4. In our model we formalize this fea- 
ture of the experiments by a switch at 7r/4; numerically 
this is accomplished by a selection criterion based on the 
value of cos 2 (ip). The model in this case encounters a 
difficulty, which is absent in a classical representation. It 
is the threshold of the PBS at the critical angle of 7r/4. 
For polarization angles close to this value the intensity in 
both channels of a PBS is equal. Theoretically, this would 
mark an event, where both polarizations, (+) and (-), 
are recorded for a single photon. Experimentally, such a 
case does not exist. Therefore we must impose a suitable 
threshold at this angle, to decide on the actual outcome. 
In our view this threshold is an experimental parame- 
ter and cannot be determined from theoretical estimates 
alone. We shall therefore use the experimental results, in 



particular the observed visibility, to determine this value. 
Formally, this reduces the range of the model PBS to the 
values cos 2 {(p) > 0.5 + As and cos 2 (<p) < 0.5 — As, where 
As is the polarizer threshold. This threshold is related to 
the intensity difference between the two channels, which 
a PBS can resolve in the experiments. A threshold of 0.1, 
for example, amounts to a resolution better than 30%. 

In all experiments of this kind decoherence of the two 
separate laser beams should play an important role. De- 
coherence in the experiments has essentially two origins: 
(i) The random deviation of the phase of laser beams 
from the phase of monochromatic waves; and (ii) the ran- 
dom deviation of the location of source and analyzers due 
to thermal vibration. In the experiments the correlations 
were independent of distance variations [l7j| . Therefore 
the laser beams are assumed to be ideal. But, as shown 
above, the angle of polarization depends on the location 
of the measurement. This means that we have to include 
thermal vibrations in the model. This is done by a ran- 
dom variation of the optical path of a photon from its 
source to the analyzer. The length of the random seg- 
ment is determined by the level of decoherence. 100% 
decoherence, for example, means that half a wavelength, 
or a phase angle of tt, is random. Correspondingly, for 
10% decoherence a phase angle of n/10 is random. 

Numerically, the procedure is as follows. 10000 pairs 
of photons are emitted from their common source. Their 
initial angle of polarization is created by a random num- 
ber generator Jl8|] , the random number [0,1] is mapped 
onto the interval [0,27r]. We add an angle of tt/2 for one 
photon of each pair. The polarization angle is set equal 
to the initial phase angle i^q. Both photons cover an op- 
tical path of X ± AX, where AX is the random segment 
due to decoherence, created by a separate random input. 
After covering the path to the analyzers, both photons 
are measured. The variables in this measurement are the 
phase angle ip (= the hidden variable) , and a polarizer an- 
gle a. We formalize the actual measurement by a switch 
cos 2 (V> -a)> 0.5 + As (+), and cos 2 (> - a) < 0.5 - As 
(-). The setting of polarizer 2 remains unchanged dur- 
ing the whole simulation. After 10000 pairs have been 
recorded, we change the angle a by 7r/100. Compared to 
the experimental practice, where the PBS remains un- 
changed, and only the speed of the optical modulator is 
altered H , the only difference should be a phase shift of 
the correlation function. In all figures we report only the 
coincidences N ++ . 

Initially we simulated an ideal measurement. The 
fields of the two laser beams are fully coherent through- 
out the distance between the two polarizers, and the ex- 
perimental devices are supposed to have ideal characteris- 
tics. The result of this simulation is shown in Fig. [l]. The 
simulation differs from the theoretical curve obtained by 
Eq. (H]). This difference is due to the digital measure- 
ments, on which our model is based. While Furry's inte- 
gral |1(| predicts a minimum of N/ 8 for the correlation 
curve, where N is the number of photon pairs, we ob- 
tain zero. The visibility (N max ~ N mm )/(N max + N min ) 
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for an ideal measurement is therefore 1, contrary to 0.5, 
obtained by the integration. The simulation also differs 
from the prediction in quantum mechanics, which is de- 
scribed by N/2 sin 2 (a - (3). The correlation curve in the 
simulation has a different functional form: while it is si- 
nusoidal in quantum mechanics, it is a saw tooth in our 
model. This difference is due to the omission of interfer- 
ence terms, or the disregard for superposition. The max- 
imum (N/2) and the minimum (zero), however, are equal 
in both cases. We also simulated the effect of the polar- 
izer threshold As on the result of ideal measurements. 
The threshold As was varied from 0.0 to 0.2. Here, we 
find that the threshold reduces the absolute yield of co- 
incidences in a simulation and increases the width of the 
minimum at the angles and 7r. This effect is equal to a 
retarded onset of the correlation function at its minimum 
position. 
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FIG. 1. Coincidence (JV++) counts of 10000 emitted pho- 
ton pairs. The ideal measurement (visibility 1.0) has the 
shape of a saw-tooth. 

The polarizers in current measurements are more than 
400 m apart [^|. Furthermore, there is no vibration 
damping or cooling to very low temperatures involved in 
such a measurement. This feature of the measurements 
is bound to cause random motion of system components. 
From surface science the range of motion without damp- 
ing can be estimated, it should be for an isolated surface 
no less than a few nanometers or more than one percent of 
the photon's wavelength. Considering that we deal with 
three coupled components and optical paths in between it 
seems safe to increase this estimate by one order of mag- 
nitude. In this case we have to include random motion 
of our system in the range of about 5-10% of the wave- 
length. This translates, in our model, into a decoherence 
rate of 10-20% (100% means that half a wavelength of the 
photon's optical path is random). Simulations with a de- 
coherence rate of 10%, 50% and 100% are shown in Fig. 
||. The interesting feature of decoherence is that it ren- 
ders the resulting distribution more sinusoidal than the 
correlation function of an ideal measurement. The fully 
decoherent simulation proves that the correlation, in our 
model, depends only on a phase connection between the 
two photons. 

If decoherence due to thermal vibrations of the system 
components is assumed to be about 10%, we may esti- 
mate the actual threshold As in the experiments from 
the experimentally obtained visibility. In the 1998 ex- 



periment the visibility was reported to be 97% || . This 
is not in keeping with an ideal combination PBS/photon 
detectors, since a threshold As = 0.0, under the con- 
dition of 10% decoherence, leads to a visibility of only 
86%. This value is also lower than the visibility of w 
reported in the 1981 experiment 17|. From the ex- 



perimental facts we have to conclude that the threshold 
As cannot be zero. This conclusion is also consistent 
with previous considerations, which established that a fi- 
nite threshold is necessary to decide on the outcome of 
a polarization measurement using a PBS. Then the only 
remaining question is the actual value of As. We only 
obtain a visibility of close to 97%, as in the 1998 experi- 
ments, if As k, 0.1. 
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FIG. 2. Coincidence counts for non-ideal measurements. 
Decoherence renders the resulting curves (10%, 50%) sinu- 
soidal, the fully decoherent curve (100%) shows that correla- 
tions in our model are only due to the connecting phase. 

In Fig. H (a) we show the result of our simulation un- 
der these conditions. It can be seen that the correlation 
curve in this case is very similar to the ideal curves from 
quantum mechanical models, it is also virtually identical 
to the curves obtained in the 1998 experiment ||. It has 
been claimed that these experiments rule out "objective 
local theories" or "realistic local theories" . The 

basis of this claim was the assertion that within both 
classes of theories the Bell inequalities are not violated 
in such an experiment. To show that our model violates 
the Bell inequalities, we simulate the counts at four se- 
lected angular positions of the polarizers a and (3 (0, 45, 
22.5, 67.5). These positions yield the maximum violation 
of Bell's inequalities in the standard framework. We per- 
formed the simulations for varying threshold values from 
0.0 to 0.2, and for a fixed decoherence rate of 10%. For 
every setting we performed 10 separate runs, each with 
10000 pairs of photons. Fig. |^ (b) gives the result of our 
simulation. The violation was computed according to the 
version of Clauser et al. (CHSH) p0|] . As shown, it in- 
creases with increasing threshold. Furthermore, the limit 
of violation is close to 2.0 (CHSH value of 3.90) in the 
final setting. For a threshold of 0.1 and a decoherence 
rate of 10%, which we found close to the experimental 
conditions, the CHSH value is 2.69 ± 0.02, as compared 
to 2.73 ± 0.02 reported in the 1998 experiments g. The 
results are therefore in excellent agreement. Finally, we 
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simulated Bell violations over the whole range of polar- 
izer angles, using the formulation employed in Aspect's 
1982 experiments Also in this simulation decoher- 

ence was 10%, the number of photon pairs was 10000. 
The result of this simulation is shown in Fig. |3| (c). It 
can be seen that it is virtually identical, for realistic con- 
ditions, to the result obtained in the experiment Q. 

Given the fact that our model is a local and realis- 
tic model of the experiments, the claim that these EPR 
experiments rule out any "realistic local theory" is only 
valid, if the experiments are in fact ideal. Which we think 
they are not. But instead of using current experimental 
results as arguments how much or how little the various 
no-go theorems apply, we shall suggest controlled experi- 
ments in future publications. There, the non-ideality can 
be related to model parameters, and these in turn allow 
extrapolating the experimental results to ideal situations 
p2fl . We think that this procedure may ultimately pro- 
vide a better understanding of the experiments than ex- 
ists at present. 
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FIG. 3. (a) Correlation curve for the most likely experi- 
mental parameters, (b) CHSH value [20] for varying threshold 
As and fixed decoherence of 10%. The Bell inequalities are 
violated in all cases where As is not zero, (c) S((j>) [21] for 
< (f> < 7r/2 and varying As. The curve for As = 0.1 is equal 
to the experimental results. 

A referee has asked whether the author actually "be- 
lieves" that this model "describes nature". I think this 



is not a scientific question. Every model, which correctly 
describes an experiment, describes nature. The only rele- 
vant question is whether one model has advantages com- 
pared to another one. Given that this model is local, it 
is realistic, and it explains the observed correlation in a 
simple manner, it has at least three advantages compared 
to the standard model in quantum mechanics. Whether 
it is correct on a more general level, depends on the re- 
sults of subsequent experimental research. 
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